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ABSTRACT 


The  mathematical  properties  of  the  following  probability 
distributions  and  their  moment  generating  functions  are 
derived: 

Weibull  distribution 
Rayleigh  distribution 
Exponential  distribution 
Normal  distribution 
Voznesensky  distribution 
Generalized  Rayleigh  distribution 

The  relative  merits  of  applying  these  distributions  to 
problems  in  ship  responses  to  the  sea,  which  is  described 
as  a stationary  stochastic  process,  are  discussed. 

In  Part  2,  the  most  promising  long-term  distributions 
derived  from  the  above  survey  are  applied  to  ship  bending 
stress  data  from  four  ships.  It  is  concluded  that  a 
numerical  solution  of  long-term  distributions,  using  either 
Weibull  or  Normal  distributions  of  the  short-term  Rayleigh 
parameters  (classified  by  weather  groups)  is  better  than 
any  explicit  function. 
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Mathematical  Properties 


1.  Introduction 

Researchers  of  the  past  decade  have  studied  the  waves  of  the  ocean’s  surface 
and  the  accompanying  wave -induced  stresses  in  a ship's  hull  by  treating  these 
quantities  as  continuous  random  functions  which  may  be  represented  by  stationary 
stochastic  processes.  The  wave  height,  in  terms  of  displacement  from  the  mean, 
may  at  any  instant  be  viewed  as  the  net  displacement  resulting  from  the  sum  of  a 
large  number  of  small  displacements,  each  of  which  is  the  result  of  a single  wave 
originating  over  a wide  area  and  which  is  independent  of  the  others.  The  central 
limit  theorem  of  probability  theory  guarantees  that  the  sum  of  a large  number  of 
independent  random  variables  may  be  treated  as  if  it  were  itself  normally  dis- 
tributed. Thus  the  process  governing  wave  height  may  be  treated  as  normally  dis- 
tributed. Since  a ship's  response  is  linearly  related  to  the  wave  height,  then 
the  response  may  also  be  viewed  as  the  sum  of  a large  number  of  independent  ran- 
dom variables  and  hence  also  treated  as  if  it  were  normal.  Thus  the  stationary 
stochastic  process  of  interest  is  assumed  to  be  normal. 

The  ship  designer  is  concerned  with  many  practical  questions  about  this 
normally  distributed  stochastic  process.  If  the  displacement  from  the  mean  is 
normally  distributed,  what  is  the  distribution  of  the  peak-to-trough  values  in 
many  reversals?  What  are  the  chances  of  exceeding  a particular  stress  level 
in  one  reversal,  or  in  an  operating  period,  or  in  a ship’s  lifetime?  What  are 
the  chances  that  the  maximum  stress  in  any  period  does  not  exceed  a certain 
level?  If  the  questions  depend  on  weather  or  sea  state  conditions,  what  are  the 
answers  to  the  above  questions  with  and  without  regard  to  weather  and  sea  state? 
What  are  the  chances  of  encountering  different  weather  and  sea  conditions?  And, 
in  short,  to  what  value  should  one  design  his  structure  with  a pre-determined  risk 
factor? 

These  questions  involve  the  distributions  of  many  random  variables  which 
are  related  to  the  normally  distributed  stochastic  process.  Numerous  data  in 
the  form  of  twenty-minute  wave  or  stress  records  have  been  laboriously  collected 
in  an  attempt  to  recognize  which  distributions  are  at  work,  and  what  are  the 
values  of  any  parameters  present.  It  is  the  purpose  of  this  part  of  the  report 
to  catalog  the  common  and  not-so-common  probability  distributions  which  re- 
searchers have  used  and  to  derive  their  significant  properties.  Hopefully,  this 
will  contribute  to  the  more  intelligent  use  of  recorded  data. 


* 

. Basic  Probability  Concepts 

Before  we  embark  on  a discussion  of  the  different  distributions,  we  shall 
briefly  discuss  some  of  the  elementary  probability  ideas  which  are  used  in  the 
comparative  discussion  which  follows. 

A continuous  random  variable  X is  described  by  two  functions  f(x)  and  F(x) 

which  are  respectively  called  the  probability  density  function  (pdf),  and  the 

cumulative  distribution  function  (cdf)  for  X.  F(x)  is  the  probability  that  the 

dF 

random  variable  X takes  on  a value  less  than  or  equal  to  x,  while  f(x)  = — . 

II  ^ 

Both  functions  are  loosely  referred  to  as  the  distribution  of  X and  f(x)  may  be 
loosely  interpreted  as  the  probability  that  X takes  on  a value  in  some  "small" 

interval  about  x.  The  probability  that  X falls  in  some  interval  [a,b]  is 

.b 

P(a  < X < b)  = / f(t)  dt 

a 

Thus,  F(x)  is  the  area  under  the  density  curve  up  to  value  x. 

1 


The  distributions  of  random  variables  can  be  characterized  by  their  mo- 
ments. The  n1^  moment  of  a random  variable  is  defined  as 

<30 

E(Xn)  - / x"  f (x)  dx, 

—CO 

and  the  first  two  moments  are  most  frequently  used  in  the  form  of  the  mean 
p = E(x)  and  the  variance  o2  = E(X2)  - E^x)2-  The  more  moments,  starting  with 
the  mean,  which  two  different  distributions  have  in  common,  the  more  similarly 
the  two  random  variables  behave. 

The  moments  of  a random  variable  are  exceedingly  useful  in  estimating  any 
unknown  parameters  from  a statistical  sample.  Suppose  one  has  a sample  x^,  X2*... 

*A  more  detailed  discussion  of  the  basic  concepts  may  be  found  in  Freund  (1), 
Feller  (2),  Meyer  (3).  See  references  listed  on  pp.  AO-41. 
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from  a population  where  the  unknown  distribution  is  thought  to  have  a certain 
form  but  with  unknown  parameters.  Then 
n 


xi  + x2 


+ x, 


may  be  taken  as  an  estimate  of  E(xn).  If  theoretical  expressions  for  E(xn)  have 

been  developed  in  terms  of  the  unknown  parameters,  then  these  estimates  of  E(xn) 

give  equations  which  may  be  solved  for  the  parameters. 

The  moment  generating  function  (MGF)  of  a random  variable  is  a frequently 

used  device  to  develop  the  moments  of  a random  variable.  If  X is  a random 

variable  with  pdf  f(x)  then  the  MGF,  M^O)  is  defined  as 

00 

M^e)  = E(ex®)  - / ex0  f(x)  dx 


Since 


X0 


l 

n=0 


dn  Mx(0) 


d0l 


tt  e 
n • 


0=0 


Mx(0) 


E(xn) 


l =T 

n=l 


E(xn) 


and  hence 


Several  other  useful  properties  of  moment  generating  functions  which  greatly  sim- 
plify the  comparison  of  random  variables  are  listed  below. 

1)  There  is  a one-to-one  relation  between  moment  generating  and  distribu- 
tions of  random  variables.  That  is  if  X and  Y are  random  variables 
whose  MGF's  My ( 0 ) and  My(0)  both  exist  and  are  equal,  then  X and  Y have 
the  same  distribution. 

2)  If  the  density  of  random  variable  X has  a parameter  A and  as  A -►  A0, 

My (0)  -y  My(0) , then  the  distribution  of  X may  be  treated  as  the  distribu- 
tion of  Y for  A close  to  A0. 

3)  If  X has  MGF  Mx(0)  and  Y - * + -*•  , then  Y has  MGF  My(0)  = ea/b0Mx(-g). 

4)  If  Xi,  X2,...,  Xn  are  independent  random  variables  with  respective  MGF's 

Mx(0) MXn(0)  and  Y » Xj_  + X2  + ...  + the  MGF  My(0)  - MXl(0) 

MX2(0) »Ve)- 

We  require  an  additional  property  of  probability  density  functions.  If  X 
and  Y are  two  random  variables  with  respective  pdf's  f(x)  and  g(y)  which  are 
related  in  a one-to-one  fashion  through  a steadily  increasing  or  decreasing  func- 
tion x = h(y)  then  1 . 1 


g(y)  “ f (h (y) ) 


dy 


3.  The  Welbull Distribution 


The  Weibull  distribution  has  been  used  of  late  by  a large  number  of  researchers 
(Nordenstrom  (4  ),  Mansour  (5  ),  Hoffman  and  Karst  (6  ) Voznessensky  (7)),  in  deal- 
ing with  engineering  problems  related  to  the  oceans  surface  and  may  be  called 
"currently  in  fashion".  The  high  flexibility  of  its  two  or  three  parameters  allow 
for  excellent  curve  fitting  to  a wide  variety  of  data,  and  many  earlier  used  distrib- 
utions may  be  made  special  cases  of  the  Weibull  distribution. 

A random  variable  is  said  to  have  a Weibull  distribution  if  its  pdf  has  the 


form 


f(x) 


a 8 x 


3-1 


The  cumulative  distribution  function  is 


F(x) 


0 


x > 0 
x < 0 


where  a > 0, 
3 > 0 


x ^ 0 
x < 0 


This  is  the  two-parameter  form  which  appears  in  most  references.  There  is  a 
three-parameter  form  of  this  distribution  which  shifts  the  left  end  of  the  density 
curve  further  leftward  and  allows  for  negative  values  of  the  random  variable  X. 
This  alternate  form  of  the  Weibull  distribution  may  be  expressed  as  density 


f(x) 


, x-a  J"1  "(¥)Y 

y ( — r > e 


and  cumulative  distribution 


F(x) 


' 1 - e 

1 0 


x ^ a 
x < a 

x ^ a 
x < a 


(Note:  Some  references  prefer  to  write  b as  V-a) 


We  shall  calculate  the  mean  and  variance  of  this  random  variable  from  its 
moment  generating  function  which  will  be  calculated  first  for  the  two-parameter 
form  and  then  converted  to  the  three-parameter  form  by  property  3 of  moment  gen- 
erating functions. 

Now  00  oo  a 

„ x0.  r x0  _ , . , r . 3-1  -axp  x0  , 

M^-C G.s  - E(e  ) e f(x)  dx  ^„J  a 3 x e e dx  , 

and  letting  y ■ ax^  we  have 

- 0(V/e 

Mx(0)  - / e “ eydy 

x 0 
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If  e 


y 1/3 

<i>  8 


is  expanded  as  £ _i 

On  • 


n/B 


we  have 


„ ...  r -n/B  0n  r -y  n/B  , 

M (6)  = l a — T Je  y dy- 

* n=0  n-  0 

00 

The  integral  / yS_1  e~y  dy  is  denoted  by  T (s)  and  is  the  well  known  gamma 
0 

function.  It  is  a generalization  of  s.'  (T( s+1)  = s.'  for  s an  interger)  and  its 
properties  are  outlined  in  the  appendix. 

Thus  the  MGF  for  the  two-parameter  form  of  the  Weibull  distribution  is 


V0)  = l fr  a n/B  r (i + ») 

n-0 


,th 


8 


and  the  nL  moments  are 


E(Xn)  = a~n/B  T(l+|) 


If  X has  the  above  3-parameter  form,  then  z ■ x-a  is  a random  variable  fitting 
the  two-parameter  Weibull  density  with  g = v and  a ~ CP—  . Thus 


-6 


n=0  ' 


Since  x = z+a,  property  3 of  moment  generating  functions  gives 
Mx(0) 


n 


e“'  l b r (1  + “) 

n=0  Y n- 


as  the  MGF  of  the  three-parameter  Weibull  distribution. 

Since  the  a in  the  three-parameter  form  is  merely  a translation,  we  shall 
restrict  our  investigation  of  properties  to  the  two -parameter  form  unless  forced 
to  use  the  three-parameter  case.  The  mean  of  the  three-parameter  form  is  merely 
the  mean  of  the  two-parameter  form  translated  by  a and  the  variances  are  the 
same. 
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First 


_1 


and 


P 

.2 


r O-  + «) 


M 1 V 

E(X2)  - E (X)2 

- 2 2 

a e [ r a + g)  - 


2 1 
r a + p ] 


Observe  that  the  parameter  a may  be  eliminated,  and  the  simple  relation 


o2  + p2 


r a + i > 

r2(i  + I ) 

6 


£ R(3) 


established  between  P and  a and  the  3 parameter.  Since  P and  o may  be  estimated 
from  a statistical  sample,  the  3 parameter  may  be  easily  estimated.  Unfortu- 
nately, the  3 may  not  be  explicitly  solved  for  and  one  must  resort  either  to  a 
numerical  solution  to  the  equation  or  to  examining  the  graph  of  R(3)  in  Figure 
1.  Note  also  that  as  3 increases,  R(3)  asymptotically  approaches  1 and  hence 
the  larger  p becomes,  the  larger  3 will  be. 

The  MGF  Mx(6)  may  be  written  as 

-1/0  n | 


Mx(9)  = I <- 


da 


or 


n: 


*^(0)  = My  (0  a-l/f 


-) 


r<-t 


> 


M -1/3  (0) 


-1/3 

where  Y is  Weibull  with  a “ 1 and  parameter  0.  Thus  X = a Y,  and  the  a appears 
to  be  playing  the  role  of  a change  of  scale  in  the  Weibull  distribution.  The  0 
parameter  seems  then  to  characterize  the  type  of  distribution  which  must  first 
be  determined  from  data.  Once  the  0 is  determined  then  the  a may  be  estimated 
either  from 

-1/0  1 

p - a r (1  + J) 


or  from  a maximum  likelihood  estimate  with  confidence  intervals  as  will  be  dis- 
cussed later. 
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The  above  discussion  implies  that  a scalar  multiple  of  a Weibull  variable 


is  still  Weibull  distributed.  This  is  S0j  for  if  Y = XX  where  X has  pdf 


f(x) 

then  Y has  pdf 


!B-1  -a  xe 
a B x e 

0 


or 


g(y) 

g(y) 


3-1 

a B ( f ) 

A 

o 

_ ' B-l  -a 

a B y e 

0 


• y x® 
-a,  f ) 


3 

y 


i 

x 


x 2l  0 
x < 0 

y j>  o 
y < 0 
y i o 

y < 0 


where  a = a/X^  . This  means  that  all  Weibull  random  variables  are  linearly  re- 
lated to  a standardized  Weibull  variable  with  a = 1.  Hence  tables  for  "stand- 
ardized variables"  at  different  g values  could  be  prepared  which  would  greatly 
simplify  use  of  this  distribution. 


4.  Maximum  Likelihood  Estimate  and  Confidence  Intervals 


« 


* 


If  X^,  X2»--*Xn  are  n independent  Weibull  random  variables  with  the  same 

fixed  8 parameter  and  identical  distributions 

...  B-l  -a  xe 

f(x)  = a 8 x e 

then  the  liklihood  function  L(a)  for  a set  of  values  xj.,  X2»...xn  is  defined 
as  the  probability  of  = x^j.^.X^  = xn.  Thus 


L(a)  = a 8 


6-1  -a  xi 


. . a 6 xn  e 


6-1  _a  xnc 


a 8 ( x^X£»  • "Xjj) 


3-1 


i=l 


The  value  of  a which  maximizes  L(a)  then  is  taken  as  the  "best"  estimate  a for 
this  sample  and  is  called  the  maximum  likelihood  estimate. 

(Note:  The  maximum  likelihood  estimate  may  or  may  not  agree  with  the  value 
obtained  through  use  of  o's  relation  to  the  mean.  These  are  two  different  es- 
timates which  reflect  two  different  criteria  of  "best"  estimate.  They  may  fre- 
quently agree  but  in  general  need  not.  The  advantage  of  the  maximum  likelihood 
estimate  in  this  case  is  that  a confidence  interval  may  be  developed  for  the 
true  a value  of  the  sampled  population  rather  than  a simple  value  with  no  error 
bounds) . 

L(a)  may  be  maximized  by  setting  ~~  = 0.  This  is  awkward,  however . Since 


d In  L(a)  „ _1 dL 

da  L(a)  da 

we  shall  solve 

d..ln  _L.fr)  , o 

da 


instead. 

In  L(a) 


n n 

« n In  a + n In  6 + (6-1)  l In  x.  - a £ x. 

i-1  i-1 


d In  L(g)  „ n _ [ 


da 


and  thus 


-1 


1 ? 
n i,  Xi 


i=l 


0 


xi 


i«l 


is  the  "best"  estimate  of  a. 

A confidence  interval  for  a,  or  equivalently  1/a  may  be  developed  as  follows. 
Consider  the  random  variable 


U = i l X* 

“ i-1  1 

where  are  independent  identically  distributed  Weibull  random  variables,  u is 
itself  a random  variable  with  a distribution  yet  to  be  determined. 

First,  if  X has  pdf 


f (x) 


.6 


Y = x has  pdf 


or 


or 


g(y) 


g(y) 


0-1  -a  X 
ax  e 


0 


g(y)  = < a 0 y 


0-1/0  -a  y 


101 


0-1/0  -ay  1-0/0 
a 0 y e i y 

P 


1 


0 

-a  y 


a e 


x >_  0 
x < 0 


y >_  0 

y < 0 

y 2.  o 

y < 0 
y >_  0 
y < 0 


Thus,  X ^is  Weibull  distributed  with  parameter  0 - 1 and  a as  in  the  sampled 
population.  Hence,  Y has  HGF 

M <8 ) - l ( £ )°  r (n+1) 

y n=0  n> 


i <:-> 

n=0 


n 


-1 


( 1 - - ) 
a 
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By  property  4)  of  moment  generating  functions,  U then  has  MGF 


«aU  <«>  ‘ d - l > 


Ma2„U  <6)  ‘ ' 29) 


by  property  3) . 

The  random  variable  W with  pdf 


2m/2  T (m/2) 


rc-2  _ w 
2 2 
w e 


w > 0 


w < 0 


is  called  chi-squared  distributed  with  m degrees  of  freedom.  It's  MGF  is 


M (0)  = (1  - 20) 

w 

Thus  we  see  that  aU  2n  is  chi-squared  distributed  with  2n  degrees  of  free- 
dom. Extensive  tables  of  chi-square  values  of  m degrees  of  freedom  are  available 

o ■ 2 

which  list  x (x,m)  for  0 _<  X _<  1 where  X (X,m)  is  such  that 


h(w)  dw  = X 


x2(a,m) 


2 

That  is  x (x,m)  is  the  value  above  which  lies  XX  of  the  area  under  the  chi- 
squared  curve.  la 


f/rrr. 


CX,rrr\) 


If  a desired  confidence  level  X were  set,  say  X « .95,  then 
P [ x2(.025,2n)  < 2naU  < x2(.975,2n)  ] - .95 


V~W  pr* 


After  rearrangement  we  obtain 


The  Rayleigh  Distribution 


A random  variable,  X,  is  said  to  be  Rayleigh  distributed  if  its  pdf  is  of 


the  form 


2x  -x2/R 

1 6 


x > 0 


x < 0 


where  R > 0. 


This  type  of  random  variable  has  been  used  by  many  researchers  (15) . Longuet- 
Higgins  ( 8 ) and  Ochi  ( 9 ) showed  that  if  the  normal  stochastic  process  govern- 
ing wave  displacement  from  the  mean  is  a narrow-band  process  (that  is  one  where 
x(t)  = A(t)  cos  (wQt  + E(t))  and  E is  small).  Then  the  envelope  A of  the  dis- 
placements in  each  small  twenty -minute  sample  of  data  is  Rayleigh  distributed 
with  parameter  R changing  from  sample  record  to  sample  record.  Since  the  peak- 
to-trough  values  of  such  a process  may  be  treated  as  twice  the  A values;  the 
values  of  the  peak-to-trough  stress  or  wave  heights  per  reversal  are  taken  as 
Rayleigh  distributed  random  variables  with 
1 v 2 

R = — l xi  , the  mean  square  of  the  sample. 
i=l 


If  we  compare  this  f(x)  to  that  of  section  3,  we  quickly  see  that  the 
Rayleigh  variable  is  nothing  more  than  the  two-parameter  Weibull  random  variable 
with  c « 1/R  and  3=2.  Thus  all  the  Rayleigh  variable's  properties  follow 


from  before. 


v = Sr  r(i  + p 


a2  = R [ r (1+1)  - r2  (1  + i)  ] 


From  the  appendix  T (3/2)  = and  T (2)  =1  so  that 


..  = - b r i -2-1 

V*  j ° Ml  ^ J • 

The  moment  generating  function  of  the  Rayleigh  variable  then  is 

= en  «/2 

Mx(e)  = I TT  R r (i  + 2) 


A ' 


This  form  of  the  MGF  is  not  too  easily  dealt  with  and  perhaps  a more  conven- 


ient form  is 

Mx(6) 


1 + 


/E  Re2/* 


6 e 


[1  + erf  (^  | ) ] 
xe- 


which  may  be  obtained  by  direct  calculation  of  E(e^). 

Here  erf  (0)  is  the  error  function  discussed  in  the  appendix. 

The  maximum  likelihood  estimate  of  R then  is 

, 11  9 

*-  i l • 

i“l  > 

which  is  the  previously  used  value.  Indeed  the  maximum  likelihood  estimate  in 

this  case  represents, being  the  sum  of  squares  of  displacements,  the  average 

energy  over  the  sample  period.  The  X confidence  interval  for  R then  is 
n 


P [ 


2 l x2 
2 i=l  xi 


R < 


2 tli 


i 


X2  (1  -X / 2 , 2n)  X2  (1  +X/2,  2n) 

This  may  then  be  used  to  develop  a confidence  for  p and  cj  through  p = /rr  R 


R[1  - it/41  if 


the  sampled  population  is  indeed  Rayleigh  rather 


n n n 

than  take  the  point  estimates  of  p % ^ I xi  and  ~ I [x^  - — I 


i=l 


i=l 


i=l 


Also,  since 

h_ 

,2 


P2  ± u2 


TU  + P )_ 

.9  . X 


p-  r* 

for  a Weibull  random  variable,  then  it  is  reasonable  to  expect  that  for  our 


(1  +P 


point  estimates  of  p and  o , we  should  have 


2 2 
q -+  Y-  * 


p2 


_.rl2)_ 

r2(B/2) 


4 

TT 


If  this  is  not  true,  then  the  Rayleigh  distribution  can  not  fit  the  data  well. 
If  it  is  roughly  true,  then  a chi-square  test  (see  Freund  (l))  will  tell  how 
accurate  the  fit  is. 
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6.  The  Exponential  Distribution 


Various  writers  have  Investigated  the  application  of  different  distributions 
to  studies  of  ship  hull  stress  statistics.  In  some  cases  they  have  found,  or 
assumed,  that  long-term  data  (several  years)  roughly  fit  an  exponential  distri- 
bution (4  ) (5).  Hence,  we  discuss  the  exponential  distribution. 


A random  variable,  X,  is  said  to  be  exponentially  distributed  if  its  pdf 
f(x) 


1 

a e n 


x > 0 

0 x <.  0 

If  we  compare  this  to  section  3,  we  see  that  the  exponential  random  variable 

is  nothing  more  than  the  two-parameter  Weibull  variable  with  g “ 1 and  a=  !/n* 
Hence  the  mean  and  variance  are 


n r (2) 


O2  = n2  { r(3)  - r2  (2)  ] 


h 


The  MGF  of  X has  been  developed  in  section  4 as 

1^(0)  * (1-pe  ) 

The  maximum  likelihood  estimate  of  n is  thus 


1 

n 


n 

I Xi 

i-1 


and  the  X confidence  interval  for  n is 

n 


l 


p t 2 i-i  X1_ 


X2  (1  - \ , 2n) 


< n < 


n 

2 I 

ial 

X2  (1  + j , 2n) 


Again,  as  in  section  5,  if  we  use  the  point  estimates  of 

1 ^ i K 2 

u ~ [x-i  + x„  + . . . + xt]  and  o2  3<  77  ( I (x.  — . 7 x.)  ] for  a 

^ k A 2 K v K K i«i 

sample  of  size  k,  it  will  only  be  reasonable  to  assume  that  an  exponential 

2  2 

random  variable  is  at  work  if  a /p  % 1 

Application  to  ship  stresses  is  discussed  in  Part  2. 
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7.  The  Normal  Distribution 


Although  the  normal  distribution  is  not  directly  applicable  to  ship  stress 
statistics,  it  will  be  shown  in  Section  2 to  be  applicable  indirectly  to  the 
problem.  Accordingly,  in  this  section  we  wish  to  list  the  properties  of 
the  normal  random  variable  and  compare  it  to  those  of  the  Weibull  random.  It 
has  been  conjectured  that  perhaps  the  two  random  are  indeed  one  and  the  same 
or  at  least  their  density  curves  are  so  close  that  they  may  be  treated  as  the 
same. 


A random  variable,  X,  is  said  to  be  normally  distributed  if  it's  pdf 
1 


f 00 


40 2 


/2tt  o 


_oo  < X < 00 


It  is  a two-parameter  distribution  whose  parameters  are  actually  the  mean  and 


standard  deviation.  If 
1 


then  the  pdf  of  z is 


-z2/2 


g(z)  - /2T-  e < Z < » 

which  is  again  a normal  distribution  which  is  referred  to  as  standard  normal. 
Although  the  cumulative  distribution  function  can  not  be  directly  calculated, 
tables  of  its  values  for  the  standard  normal  case  are  readily  available,  and 
these  may  quickly  be  related  to  any  normal  random  variable. 

We  first  obtain  the  MGF  of  the  standard  normal  variable  z. 


M (6)  - 


\ f z6  -z2/2  . 

■ I eedz 


[ -1/2  tz2-2z6] 


dz 


We  complete  the  square  of  the  exponent  and  have 


M (e)  - e 


e2/2 


1 


-1/2  [z2-2z0+02] 


dz 


Since  the  general  normal  random  variable  X may  be  written  as  X = o(r+  — ), 


property  3 again  yields 
M^(6)  - 


6^2 

2 


U0 


The  normal  distribution  is  the  most  frequently  used  and  most  important 
probability  distribution  because  of  the  central  limit  theorem  which  guarantees 
that  the  sum  of  a large  number  of  independent  random  variables  may  be  treated 
as  if  the  sum  were  normally  distributed.  Property  4 of  moment  generating  func- 
tions enables  us  to  apply  the  central  limit  theorem  to  approximate  many  common 

random  variables  via  the  normal  distribution.  For  example  the  aforementioned 

-m/2 

chi-squared  distribution  with  m degrees  of  freedom  had  MGF  of  (1  - 20) 

When  m=2n,  this  is  the  product  of  n MGFs  of  (1  - 20) which  is  the  MGF  of 
an  exponential  random  variable.  Thus,  the  chi-square  variable  may  be  viewed 
as  the  sum  of  n independent  exponentially  distributed  random  variables  and 
hence,  via  the  central  limit  theorem,  approximated  by  normal  random  variable 
with  the  same  mean  and  variance. 

In  order  to  see  how  well  the  normal  and  Weibull  distributions  agree,  we 
must  first  decide  which  form  of  the  Weibull  to  use.  Since  the  normal  may  take 
on  negative  values,  the  three-parameter  form  initially  seems  more  appropriate. 
However,  since  the  discrepancy  between  the  two  forms  applies  to  energy  states 
which  are  always  positive,  the  two-parameter  form  seems  to  be  the  one  to  com- 
pare to  the  normal.  Thus  we  compare  the  moment  generating  functions: 


and 


2 2 

0 o + we 

e 2 


l 

n**0 


n 

6 


nj 


-n/8 

a 


r(i  + f> 


(normal) 

(two-parameter  Weibull) 


n 

Even  when  the  normal  MGF  is  expanded  as  a series,  the  presence  of  T(1  + *g)  in 
the  Weibull  MGF  make  the  two  functions  sufficiently  different  that  it  appears 
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there  is  no  hope  of  making  them  agree.  Thus,  the  normal  and  two-parameter 
Weibull  distributions  are  two  distinctly  different  distributions! 

Let  us  still  examine  the  graphs  of  these  two  density  functions  over  the 
range  of  values  usually  dealt  with.  The  mean  and  variances  of  a probability 
density  function  characterize  the  centering  and  spread  from  the  center  of  the 
probability  weight.  The  higher  moments  (3rd,  4th...  etc.)  of  the  distribution 
characterize  the  skewness  of  the  density  about  the  mean.  (That  is  the  tendency 
to  flatness  in  the  density  curve  on  one  side  of  the  mean.)  If  two  probability 
density  curves  are  to  be  roughly  the  same  then  it  is  clear  that  in  addition  to 
their  means  and  variances  being  equal,  as  many  higher  moments  as  possible  must 
be  also  equal.  Thus,  if  the  normal  and  Weibull  are  approximately  equal,  we  must 
have 


a 1/e  r (1  + £ ) 


2,2 
a + p 


.-2/  sra+i) 


(1st  moment) 
(2nd  moment) 


3po2  + p3  « a r(l  + “■  ) (3rd  moment) 

etc. 

When  the  first  two  equations  are  inserted  into  the  third,  we  get  the 
condition 

r(i+“)  [ 3 r (l  + - ) - 2 r2  (l  + ) ) = rU+f> 

or  o i o ^ 

o = H(6)  = rd  + f ) - rd  + g ) [ 3 rd  + g ) - 2 r2(i  + g ) ] 

A value  of  B which  makes  H(B)  ^ 0,  in  addition  to  satisfying  the  first  two 
equations,  will  produce  a Weibull  density  quite  similar  to  a normal  density. 

An  investigation  of  H(($)  (see  Figure  2)  shows  that  H(g)  * 0 for  B^.S. 
However,  for  3 B 10,  |h($)|<5  x 10  . Thus,  it  appears  that  a Weibull 

distribution  will  behave  "like"  a normal  distribution  for  3 <_  B <.  10,  and  that 
the  closest  approximation  will  be  when  g ^ 3.5.  It  is  reasonable  to  ask  how 
the  4th  and  higher  moments  behave  when  H(g)  ■ 0.  Unfortunately,  when  H(g)  * 0, 
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the  Ath  moments  do  not  agree.  Hence,  the  Weibuil  will  not  tend  to  a normal 
exactly,  but  it  can  be  made  to  fit  up  to  the  3rd  moment.  The  accompanying 
graphs  (Fig.  3)  illustrate  how  good  this  fit  is. 


8.  The  Voznesensky  Distribution 

A brief  summary  of  the  Voznesensky  distribution  may  be  found  in  Ochi  (9). 

/ 

This  distribution  was  developed  from  an  empirical  standpoint  using  the  random 
samples  of  peak  values  instead  of  a spectral  analysis. 

Voznesensky  defined  his  distribution  through  its  moments.  First,  he  re- 
quires a random  variable,  X,  where  k^  moments  were  E(X^)  = r ( 1 + ^ ) for 
some  still  to  be  determined  n.  Next  the  random  variable  governing  the  distribu- 
tion of  peak  values  per  reversal  was  taken  as  Z - X/T(l  + . The  n could 

then  be  determined  using  point  estimates  of  \i  and  o from  the  equation 

1 = rq  + 2/n)  - r2q  + i/n) 

M r (1  + l/n) 

If  we  compare  the  moments  with  those  of  the  Weibuil  variable  in  section  3, 
we  see  that  Voznesensky’s  X variable  is  merely  a two-parameter  Weibuil  variable 
with  a * 1 and  6 = n.  His  condition  for  determining  n is  the  condition 

q-  ± -V2-  = R(  8) 

P2 

which  was  discussed  before.  Hence  we  shall  not  further  discuss  this  special 
case  of  the  Weibuil  distribution. 
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9.  The  Generalized  Rayleigh  Distribution 


The  earlier  work  of  Rice  (10),  is  summarized  in  Ochi  (9).  Rice  examined 
the  stochastic  process  x(t)  governing  the  displacement  from  the  mean.  Assuming 
x(t)  is  normal,  the  distributions  and  correlations  of  x(.t)  and  x(t)  were  de- 
rived. Let  f(x(t),  x(t),  x(t))  be  the  joint  density  of  these  three  random 
variables.  Maxima  of  x(t)  occur  when  x(t)  is  zero  while  x(t)  is  negative,  and 
probability  of  peaks  of  x(t)  exceeding  level  £ may  be  obtained  by  integrating 
the  joint  density  function  with  x(t)  = 0,  over  the  range  £ < x(t)  and  x(t)  < 0. 
The  chief  advantage  of  this  approach  is  the  appearance  of  the  spectral  band- 
width parameter  e in  the  distribution. 

When  the  above  process  is  finished,  we  arrive  at  the  following  density 
function  for  X,  the  maxima  of  the  stochastic  process. 


f(x) 


1 

/2tt 


e 


-x2/2e2 


+ A - e2 


-x2/2 

xe 


j € ^ 

.00 


for 

-°°  <x<°°  , 0<e<]_ 

A random  variable  with  this  density  is  said  to  have  the  generalized  Rayleigh 
distribution.  The  name  follows  from  the  fact  that  as  £•+■  0,  and  hence  the 
process  becomes  narrow  band,  the  distribution  approaches  that  of  a Rayleigh  ran- 
dom variable  with  R “ 2.  (Note:  In  the  above  density  function  we  have  already 
divided  X by  its  rras  value,  and  hence  the  additional  parameter  does  not  appear 
in  the  density  function  and  R is  constant).  As  e-*-l,  however,  and  the  band  width 
widens,  the  density  function  approaches  the  standard  normal  case.  A graph  of 
f(x)  for  various  t values  is  in  Figure  4. 

Ideally,  we  would  like  the  MGF  for  the  random  variable.  Although  the  form 
of  the  density  is  quite  involved  and  initially  seems  to  prohibit  a simple  in- 
tegral evaluation,  several  pages  of  integral  transformations  listed  in  the 
appendix  yield 
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Mx(6) 


|2e2/ 2 + &A  - e2  /V2  e02/2  [ 1 + erf  ( -BiGL  ) J, 


As  expected,  we  see  that  as  e-K),  the  MGF  becomes 

2 


Mx(6) 


1 + d/ir/2  e0  /2  [ 1 + erf  ( 9 )] 

/F 


which  is  that  of  a Rayleigh  random  variable  with  R=2.  As  e->T,  the  MGF  becomes 


Mx(e) 


,6 


2/2 


which  is  that  of  the  standard  normal  random  variable. 

If  the  MGF  Mx(0)  is  expanded  as  a powers  series  in  0,  we  have 


Mx(9) 


1 +/tt7"2  A - A2  0 + [2  - e2]  0 2/2 ! + 3/it/2  - e2  03/3.’  +. 


from  which  we  may  easily  read  off  E(X)  and  E(XZ).  Hence, 
p = - e2  /ir/2 

a2  = 1 + (1  - e2)  (1  - tt/2). 

This  of  course  agrees  with  the  mean  and  variance  of  the  two  limiting  cases. 

Before  we  compare  this  distribution  to  the  Weibull  case,  we  shall  make 
some  interesting  observations.  The  MGF  M^  (0)  above,  may  be  written  as 


Mx(0) 


m7  ( -- 


/2 


where  Z is  a random  variable  with  MGF 


Mz(9)  = 

Hence,  by  property  3 


02e2/d-e2) 


— -0 2 
1 + /it  0 e [1 


+ erf  (0)  ] | 


'1-t 


Z. 


/2 


But  M7(0)  is  the  product  of 


and 


Mw(0) 


02e2/(l-e2) 


My(e)  - 1 + /it  0 e°  [1  + erf  (0)]. 
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•w-i 


~V~w  - 


M^(0)  is  the  MGF  of  a normal  random  variable  with  0 and  a2  = 2e2/(l-€2). 
Hence  W has  density 


h(v) 


S\  - £ 2 

S%7TT 


_ w2(l-e2) 
A£  2 


_ oo  < ^ < oo 


My (9)  is  clearly  the  MGF  of  a Rayleigh  random  variable  with  R=2  and  hence 


density 


gty)  = 


y e 
0 


-y  / 2 


y >_  0 

y < o 


Thus  by  property  4,  the  random  variable  X may  be  written  as  the  sum  of 
two  independent  variables  U and  V 


9 2 

One  variable,  U,  is  normal  with  p = 0 and  o = e . The  other,  V,  is 
Rayleigh  with  p =/ir  /2  /l-f.  2 and  o2  = (2  - it  / 2 ) ( 1 - e2).  The  Chebyshev  Theorem 
states  that  ^ny  random  variable  has  no  more  than  a 1/n  change  of  falling  more 
than  n cr  distance  from  its  mean.  Thus,  for  e near  zero,  the  U values  have  little 
effect  on  X,  and  for  e near  1,  the  V values  have  little  effect.  Thus,  the  more 
narrow-band  the  stochastic  process  becomes,  the  more  the  distribution  of  the 
maxima  becomes  Rayleigh. 

If  we  compare  the  generalized  Rayleigh  case  to  the  Weibull  case,  we  immed- 
iately see  that  the  three-parameter  Weibull  distribution  must  be  used  since 
X may  take  on  negative  values.  Since  the  form  of  MGFs  are  again  quite  distinct 
as  in  section  7,  there  is  no  hope  of  showing  that  the  generalized  Rayleigh  var- 
able  is  really  a three-parameter  Weibull  variable  withy,  b,  a related  to  e. 
However,  since  the  Weibull  distribution  does  have  three  parameters,  y,  b,  a, 
these  may  be  selected  as  functions  of  e to  require  that  the  first  three  moments 
of  the  two  distributions  fit.  This  will  make  the  mean,  variance  and  skewness 
coincide  again,  and  hence  the  two  curves  would  be  quite  similar  to  the  eye. 
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Thus  we  want 

/tt/2  A -e2  = b r(l  + ) + a = E(X) 

2 - e2  = a2  + 2ab  T(1  + ^ ) + b2  T(1  + ^ ) *=  E(X2) 

3/172  X - e2  = b3  T(1  + ^ ) + 3 ab2  T(1  + | ) + 3 a2b  T(1  + 

+ a3  = E(X3) 

Now  if  b = /l  - e2  B , a =/l  e2  a 
The  equations  become 

/7/2  = B T (1  + ^ ) + a 

1 + _i — = a2  + 2aB  T(1  + I ) + B2  T (1  + - ) 

l-e2  y Y 

3/72  ■=  B3  T(1  + -^  ) + 3 a B2  T(1  + -^  ) + 3 a2  B T(1  + 1 ) + a3 

These  simplify  to 


Thus  a "good  fit"  of  the  Weibull  density  to  the  generalized  Rayleigh  distribu- 
tion may  be  achieved  with  the  dependence  on  e,  the  spectral  width,  clearly 
brought  out. 
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10. 


Further  Welbull  Properties  (two-parameter) 


A)  The  Mode: 

The  mode  of  a random  variable  X Is  that  value  which  is  most  likely  to  occur 
and  will  be  the  most  frequently  observed  value  in  a large  sample.  The  mode  will 
occur  when  the  pdf  f(x),  or  equivalent.  In  f (x)  achieves  its  maximum.  Now 

In  f(x)  = In  a B + (B  - 1)  In  x - a xB 


d In  f(x) 
dx 

which  vanishes  when 


(0  - 1) 


i - ctB_  xe_1  i 
x B-l  1 


x 


1/8 

] 


B)  The  p-th  value  Xpi 

Xp  is  the  value  above  which  lies  p%  of  the  probability  weight,  or  above 
which  p%  of  the  values  in  a large  sample  will  lie.  Thus 


00 

Jp  £«> 

dt  * p 

or 

B 

-aXp 

e « 

P 

Hence 

B 

axp 

In  (1/p) 

and 

XP  “ 1 

In  q/e)  ] 

a 

C)  Order 

Statistics : 

If  a sample  of  n independent  random  variables  is  taken  a large  number  of 
times,  a histogram  of  the  maximum  values  can  be  developed.  The  maximum  value 
of  a sample  of  sire  n is  itself  a random  variable,  Xmax,  with  a distribution. 
Now  requiring  X^^  < X is  equivalent  to  requiring  that  each  of  the  n sampled 
random  variables  XA  <_  x for  i - l,2,...,n.  Thus  if  F(x)  is  the  edf  of  the 
sampled  population,  then 
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a 


P<Xmaxix>  “ F(x) 
is  the  edf  of  Xmax  and 
n-1 

n F (x)  f(x) 
is  the  pdf. 

Thus  for  the  Weibull  random  variable  the  density  function  of  \tax  is 


r -ax^  . g-1  -ax® 

a g n I 1 - e ] xp  e 


0 < x 


This  can  be  maximized  easily  to  find  its  mode.  This  will  be  the  most  frequently 
occurring  maximum  value  observed  in  many  repeats  of  n samples.  The  probability 
of  a maximum  exceeding  this  value  is  so  large  as  to  be  of  little  design  interest. 
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11.  Numerical  Solution 

Finally,  mention  should  be  made  of  the  possibility  of  obtaining  a long- 
term distribution  by  numerically  summing  a large  number  of  short-term  distri- 
butions, such  as  the  Rayleigh  (6).  In  this  case  it  is  unnecessary  that  the 
resulting  distribution  fit  any  particular  function.  Of  course,  the  method 
can  be  irefined  by  using  the  generalized  Rayleigh  distribution  for  the  short- 
term data,  having  an  additional  bandwidth  parameter,  e,  as  discussed  in 
section  9 ( 6 ) • 

Available  ship  stress  data  are  usually  in  the  form  of  short-term  (20  to 
30  minutes)  records  taken  automatically  ever  four  hours.  If  the  peak-to-trough 
stress  variations  are  assumed  to  be  Rayleigh  distributed,  then  the  individual 
records  can  be  characterized  by  their  Rayleigh  parameters,  R (rms  values). 
Hence,  one  way  to  proceed  is  to  determine  a suitable  distribution  function 
to  describe  the  R,  which  may  be  treated  as  a random  variable  with  pdf  g(R). 

If  X is  the  measured  stress  or  wave  height,  the  Rayleigh  distribution  obtained 
from  each  record  may  be  considered  a conditional  distribution  of  X given  R 
which  we  denote  by  f(x/R).  Then  f(x/R)g(R)  is  the  joint  density  of  X and  R, 
and  the  density  of  X alone  regardless  of  R may  be  obtaineu  from: 

r *(*/*>  8(R)  dR 
0 

This  is  what  is  referred  to  as  the  long-term  distribution. 

Examples  will  be  presented  in  Part  2. 

4 

The  long-term  distribution  of  stress  or  wave  heights  per  reversal  has 
also  been  obtained  as  follows  (11),  (12),  (13),  (14).  First  the  data  are 
classified  according  to  weather  severity.  Within  each  weather  group  the 
twenty-minute  records  each  yield  a Rayleigh  distribution  but  with  different 
R values.  The  R values  are  then  taken  as  a sample  of  the  distribution  of  ship 
stresses  within  that  weather  group,  and  long-term  distributions  are  obtained 
for  each  weather  group.  The  results  for  all  weather  groups  are  then  weighted 
according  to  chances  of  weather  occurences  and  added.  Examples  are  given  in 
Part  2. 
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This  approach  has  the  advantage  of  relating  the  observed  data  to  the 
physical  cause — the  sea  conditions — rather  than  relying  on  the  adoption  of 
a particular  distribution  function  that  happens  to  match  the  data  at  low  N 
values.  This  is  claimed  not  only  to  result  in  more  reliable  extrapolation 
of  the  data  to  large  values  of  N but  permits  comparison  of  ships  on  different 
services  by  reducing  results  to  the  same  "standard"  or  typical  weather 
conditions. 

Within  a given  weather  condition,  Nordenstrom  has  claimed  (4)  the  parameter 

of  short  term  Rayleigh  distribution  follows  a Weibull  distribution,  while  Webb 

researchers  (11) , (6)  have  claimed  that  the  normal  distribution  is  at  work. 

If  one  examines  a sample  of  the  parameter  X^,  using  the  point  estimates 

2 

of  p and  a to  estimate  a and  6 discussed  above,  the  g either  does  or  does  not 
fall  in  the  range  where  the  H(g)  = 0 and  where  the  two  distributions  are 
close.  The  R(g)  in  section  3 which  determines  6 must  be  close  to  R(3.5)  * 1.1002, 
if  a normal  is  also  to  fit  the  data  fairly  well.  Since  3 £ g <_  10  seems  to 
give  fairly  normal  approximations  and  (see  graph)  since  in  this  range  1.132 
R(g)  <_  1.014,  it  appears  that  a normal  and  Weibull  distribution  may  be  quite 
similar  over  a broad  range  of  a. 

A few  examples  run  on  the  weather  data  in  the  appendix  of  SSC-196  seem  to 
indicate  that  the  weather  groups  III,  IV,  V each  yield  a g within  this  range, 
but  that  the  data  for  groups  I,  II  yield  g falling  slightly  outside  the  broad 
range.  The  lumped  weather  data  also  does  not  fall  within  the  g range.  Since 

R(g)  = o2  + p2  , the  "best"  g will  occur  when  o2/p2  *v  .1002,  it  seems  reason- 
P2 

able  that  the  higher  weather  groups  with  the  larger  p values  will  yield  g's 
falling  within  a range  where  the  normal  and  Weibull  densities  are  close. 
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12.  Comments  and  Conclusions 


The  application  of  the  Weibull  distribution  to  problems  in  ship  research 
appears  particularly  promising.  Under  the  assumptions  that  a normal  "stochastic" 
process  is  initially  at  work,  a reasonable  assumption  in  light  of  the  central 
limit  theorem,  the  distributions  of  the  variables  of  interest  may  very  well  be 
special  cases  of  the  Weibull  distribution.  Those  which  are  not  directly  Weibull 

may  be  such  that  their  density  curves  will  be  fit  quite  well  by  the  Weibull 
density,  sd  that  questions  based  on  areas  under  these  curves  may  be  approxi- 
mately the  same.  Accordingly,  Part  2 of  this  report  will  examine  actual  full- 
scale  statistical  data  from  ship  research  to  determine  the  applicability  of  the 
Weibull  distribution,  as  well  as  other  approaches  discussed  here. 


Part  2 


Application  of  Long-Term  Distributions 


1.  Introduction 

It  is  the  objective  of  this  part  of  the  report  to  show  the  application 
of  the  most  promising  forms  of  long-term  distribution  presented  in  Part  1 
to  a marine  system  response  to  stochastic  inputs.  Specifically,  ship  hull 
stresses  in  irregular  waves  at  sea  represent  a phenomenon  of  great  interest 
to  ship  designers  for  which  considerable  full-scale  statistical  data  are 
available. 

A review  of  published  data  revealed  that  histograms — and/or  cumulative 
statistics — were  available  for  the  following  ships: 

Esso  Malaysia,  190,000-ton  deadweight  tanker 

R.G.  Follis,  66,500-ton  deadweight  tanker 

Fotini  L,  74,000-ton  deadweight  bulk  carrier 

Wolverine  State,  15,000-ton  general  cargo  ship 

The  data  describe  the  total  population  of  cycles  of  peak-to-trough  midship 
bending  stress  obtained  over  periods  of  two  to  three  years.  They  were 
recorded  automatically  for  periods  of  20-30  minutes  every  4 hours. 

It  was  hoped  to  obtain  similar  data  for  the  SL-7  container  ship 
Sealand  Maclean,  but  it  was  found  that  the  available  data  were  not  in 
sufficient  detail. 

In  the  case  of  the  first  ship  mentioned  above,  the  tanker  Esso 
Malaysia  computer  cards  were  available  at  Webb  Institute  for  the  entire  \\ 
years  of  data  collection.  (These  were  obtained  from  Teledyne  Engineering 
Services  in  connection  with  a project  at  Webb  sponsored  by  the  American 
Bureau  of  Shipping  (15).)  Hence,  it  was  decided  to  use  these  cards  to 
make  a completely  independent  statistical  analysis  of  the  data,  rather 
than  to  rely  on  the  previously  published  results.  The  results  of  this 
analysis,  presented  in  the  next  section,  were  found  to  agree  quite  well 
with  the  published  figures  (15). 
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It  should  be  noted  that  for  the  analysis  of  ship  stresses  the  longest 
time  for  which  data  have  been  collected  on  any  one  ship  is  about  three 
years,  whereas  a typical  ship's  lifetime  is  20  to  25  years.  Hence,  the  real 
problem  is  to  obtain  a cumulative  distribution — or  probability  model — that 
can  be  extrapolated  by  a time  factor  of  10  or  more.  Basically,  there  are 
two  approaches,  one  to  find  a long-term  distribution  that  describes  all  the 
data  (stress  reversals)  and  the  other  is  to  assume  that  all  short-term  data 
(20-minute  records)  can  be  described  by  Rayleigh  distributions  and  that  the 
long-term  distribution  can  be  obtained  by  summing  up  these  short-term 
distributions.  In  the  latter  case  it  is  not  necessary  to  obtain  a formal 
long-term  distribution  function. 

In  addition  to  the  Esso  Malaysia  analysis  presented  in  the  next  section, 
long-term  trends  for  the  other  ships  listed  above  will  be  presented  in  the 
following  section.  In  all  cases  the  applicability  of  the  Weibull  distribution 
discussed  in  Part  1 will  be  considered. 
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2.  Esso  Malaysia  Analysis 


Stress  Data 

Bending  stress  data  from  24  instrumented  voyages  of  the  supertanker  Esso 
Malaysia  were  analyzed.  Stresses  were  recorded  in  20-minute  records  ("intervals"), 
normally  one  interval  per  4-hour  watch.  The  reduced  data  population  for  the 
24  voyages,  after  discarding  intervals  containing  unusable  data,  consisted  of 
3589  analyzable  records.  Information  recorded  in  each  interval  included: 

1.  Ship,  wind  and  wave  direction  information. 

2.  Counts  of  number  of  stress  reversals  (peak-to-peak)  falling 
within  each  of  16  stress  ranges. 

3.  Total  number  of  stress  reversals  in  the  interval. 

4.  RMS  peak-to-peak  stress  for  the  interval. 

5.  Maximum  single  peak-to-peak  stress  value  in  the  interval. 

Our  analysis  of  these  data  has  been  concentrated  mostly  on  the  distribution 
of  the  stress  reversals,  item  2 above. 

A total  of  785,511  stress  reversals  were  counted  and  grouped  into  16 
stress  ranges  within  each  record  by  a Probability  Analyzer,  a small  computer 
used  by  Teledyne  Materials  Research  to  generate  the  statistical  characteristics 
(digital)  of  the  recorded  intervals  (analog).  Results  were  recorded  on  3589 
computer  cards  (one  for  each  interval) , and  it  was  this  set  of  cards  that  was 
further  processed  at  Webb  Institute  for  this  project. 

The  stress  ranges  into  which  the  Probability  Analyzer  groups  the  digitized 
stress  reversals  could  be  adjusted  to  various  levels  depending  on  the  magnitude 
of  the  maximum  measured  stresses.  In  the  Esso  Malaysia  analysis,  most  tape  reels 
had  maximum  stresses  below  8000  PSI,  so  the  analyzer  was  adjusted  to  16  ranges 
of  500  PSI  each.  For  reels  containing  higher  stresses,  ranges  of  750  PSI  each 
(maximum  stress  12,000  PSI)  or  1000  PSI  each  (maximum  stress  16,000  PSI)  were 
used.  The  flexibility  enabled  the  stress  ranges  to  be  chosen  to  suit  the  data 
without  placing  excessive  numbers  of  reversals  in  any  one  range.  The  grouping 
of  data  from  the  Esso  Malaysia  is  summarized  in  Table  I. 
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Table  I 

Summary  of  Stress  Analysis  — Esso  Malaysia 


No.  Reels  No.  Intervals 

48  2727 

10  648 

3 214 

3589 

To  get  the  whole  population  into  one  histogram,  the  numbers  of  stress 
reversals  falling  within  each  stress  range  had  to  be  summed.  This  required 
regrouping  into  constant  stress  ranges.  This  was  done,  approximately,  by 
plotting  cumulative  curves  of  the  three  groups  of  data  and  adding  them, 
as  in  Figure  5.  The  resulting  cumulative  curve  was  then  read  back  at  stress 
level  ranges  of  500  PSI,  and  the  cumulative  number  of  reversals  converted 
into  "percent  equal  to  or  less  than",  by  dividing  number  of  occurrences  by 
the  total  number  of  reversals,  N = 785,511.  Division  by  N + 1 instead  of  N 
yielded  a "plotting  position"  for  the  probability  paper  analysis  which  makes 
the  first  (lowest)  and  last  (highest)  measurements  symmetrical  with  respect 
to  the  0%  and  100%  levels,  respectively.  This  adjustment  enabled  the  most 
extreme  data  point  to  be  plotted  (4),  The  resulting  cumulative  distribution 
and  plotting  positions  are  tabulated  in  Table  II. 

Distribution  of  the  Stress  Reversals 

The  probability  of  exceeding  different  levels  of  stress  actually  measured 
on  the  ship  can  be  determined  from  Table  II.  It  is  the  "percent  greater  than" 
the  given  stress  level,  that  is: 

Q(X  > Xj)  = 1 - P(X  < Xj) 

where 


Q(X  > Xj)  * probability  of  exceeding  Xj 
P(X  < X<)  - probability  of  not  exceeding  X-t 


Probability  Analyzer  Stress  Group 

Max.  Set  ting,  PSI  Range,  PSI 


8000 

500 

12  000 

750 

16000 

1000 

" ^ 


Table  II 

Cumulative  Stress  Reversal  Histogram,  Esso  Malaysia 


Stress  Level, 
PSI 

Cum.  No. 
Reversals 
= or  < 

Cumulative 

Percent 

Prob.  Paper 
Plotting 
Position 
Percent 

Probability 

of 

Exceedance 

-lo®10P 

0 

0 

0 

0 

500 

306,700 

39.04 

39.04 

0.6096 

0.21 

1,000 

564,600 

71.88 

71.88 

.2812 

0.55 

1,500 

680,740 

86.66 

86.66 

.1334 

0.87 

2,000 

733,900 

93.43 

93.43 

.0657 

1.18 

2,500 

757,600 

96.45 

96.45 

.0355 

1.45 

3,000 

769,780* 

98.00 

98.00 

.0200 

1.70 

3,500 

776,300 

98.83 

98.83 

.0117 

1.93 

4,000 

780,900 

99.41 

99.41 

.0059 

2.23 

4,500 

782,400 

99.60 

99.60 

.0046 

2.40 

5,000 

783,400 

99.73 

i 99.73 

.0027 

2.57 

1 

5,500 

784,420 

99.86 

99.86 

.0014 

2.85 

6,000 

784,876* 

99.919 

99.919 

.00081 

! 3.09 

6,500 

785,070 

99.944 

99.94 

.00056 

1 3.25 

7,000 

785,295 

99.973 

99.972 

.00027 

! 3.57 

7,500 

785,369 

99.982 

99.982 

.00018 

3.74 

8,000 

785,411 

99.987 

99.987 

.00013 

3.89 

9,000 

785,455* 

99.9929 

99.9927 

.000071 

4.15 

10,000 

785,477 

99.9957 

99.9955 

.000043 

4.37 

11,000 

785,492 

99.9976 

99.9975 

.000024 

4.62 

12,000 

785,506* 

99.9994 

99.9992 

.000006 

5.20 

13,000 

785,510* 

99.9999 

99.9997 

.000001 

5.90 

** 

* 

13,800 

785,511 

100. 

99.9999 

I 

* 

These  are  exact.  Others  are  close  approximations,  not  exact  counts,  because  of 
the  varying  stress  level  ranges  used  in  the  original  analysis  (see  text). 

** 

Highest  single  value  of  stress  was  13,800  PSI. 


-33- 


The  probability  of  exceedance  thus  derived  from  the  stress  reversal  histogram 
is  plotted  in  Figure  6 . For  a long-term  probability  of  exceedance  of  stress 
levels  beyond  those  actually  measured,  the  curve  in  Figure  6 must  be  extra- 
polated to  higher  stress  levels  and  smaller  probabilities.  To  do  this,  a 
probability  distribution  model  must  be  determined  which  fits  closely  the 
observed  data  and  therefore  gives  a reliable  estimate  of  predicted  extreme 
values. 

Most  attempts  at  fitting  a distribution  to  the  entire  population  of 
stress  reversals  have  not  produced  satisfactory  fits  at  extreme  stress  levels. 

However,  it  appears  from  Part  1 that  one  type  of  extreme  value  distribution, 
the  Weibull  distribution,  may  satisfactorily  describe  the  entire  population 
of  stress  reversals.  This  hypothesis  was  tested  for  the  Esso  Malaysia  data 
by  plotting  the  cumulative  stress  percentages  on  Weibull  probability  paper, 
on  which  a true  distribution  plots  as  a straight  line.  Figure  7 shows  such 
a plot,  using  the  stress  levels  and  plotting  positions  given  in  Table  II. 

It  is  clear  from  this  plot  that  the  entire  population  of  stress  reversal 
data  is  not  well  fited  by  a Weibull  distribution,  since  the  plotted  points 
describe  a curve  rather  than  a straight  line.  The  question  remains  whether 
an  approximate  straight  line  fit  can  be  determined  such  that  the  corresponding 
Weibull  distribution  will  define  a satisfactory  long-term  curve  to  extrapolate 
the  measured  data.  Two  visual  fits  were  tried,  as  shown  on  the  figure.  Line 
1 was  fitted  to  the  entire  range  of  data,  and  Line  2 favored  the  high- stress 
data,  where  the  trend  was  nearly  linear.  In  both  cases  the  Weibull  parameters 
were  estimated  from  the  plotted  straight  lines.  The  cumulative  Weibull 
distribution  is: 

-ax& 

Q(x)  = e = probability  of  exceedance,  where 

x = stress  level 

4 
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The  Weibull  paper  is  constructed  with  linear  scales  of  In  x as  abscissa 
(horizontal  scale),  and  In  In  ^ as  ordinate  (vertical  scale).  A straight 
line  on  the  Weibull  paper  has  slope  = 3 and  the  intercept  with  the  zero 
axis  of  In  x is  a.  That  is,  the  straight  line  has  the  equation 


in  In  ^ 


6 In  x + In  a 


where  a and  3 are  the  Weibull  parameters. 


The  parameters  determined  from  the  Weibull  fits  are  as  follows: 


Line  1 a = 0.00060  3 = 1.091 

Line  2 a = 0.00512  3 = 0.826 


Each  of  these  lines,  in  addition  to  the  actual  curve  from  the  histogram,  was 
used  to  predict  long-term  trends.  The  long-term  curves  are  plotted  in  Figure 
8,  which  shows,  as  expected,  a poor  fit  for  Line  1.  Line  2 is  a better  fit 
at  high  stress  levels,  but  since  the  number  of  stress  reversals  at  these  levels 
is  so  small,  the  reliability  of  the  Weibull  curve  as  an  extrapolator  is 
questionable. 

A third  Weibull  fit  was  tried,  using  calculated  estimates  of  a and  3 
determined  as  explained  in  Part  1,  rather  than  from  a Weibull  paper  plot. 

The  resulting  long-term  curve  (a  = .00041,  3 ■ 1.192)  was  found 

to  be  entirely  inadequate  to  represent  the  data.  This  numerical  technique 

cannot  "bias"  the  fitted  line  toward  any  particular  part  of  the  histogram. 

All  three  fitted  Weibull  lines  were  then  subjected  to  chi-squared  goodness 
of  fit  tests  to  the  measured  data.  In  all  cases  the  results  showed  poor 
to  very  poor  correlation  with  the  data.  It  is  concluded  that  the  Weibull 
distribution  is  not  an  acceptable  model  of  unstratified  stress-reversal  data. 


3.  Analysis  of  Data  for  Other  Ships 

The  available  stress  reversal  measurements  of  several  other  ships  were 
also  tested  for  Weibull  fits,  and  the  results  are  shown  in  Figure  9.  The 
ships  are  the  R.G.  Follis,  the  Fotini  L.  and  the  Wolverine  State.  It  can 
be  seen  that  the  Follis  and  Fotini  stresses  exhibit  the  same  characteristic 
curve  as  that  found  for  the  Esso  Malaysia  data  when  plotted  on  Weibull  paper, 
and  the  fit  to  a straight  line  is  just  as  unsatisfactory.  Only  the  Wolverine 
State  data  are  well  described  by  a Weibull  distribution.  Why  this  should 
be  so  is  not  known,  but  it  is  clear  that  the  Weibull  distribution  cannot 
be  assumed  to  describe  bending  stress  distribution  for  all  types  of  ships. 

) 

Distribution  of  Rayleigh  Parameters 

Since  the  stress  reversals  do  not  seem  to  be  well  described  by  a Weibull 

distribution,  it  is  of  interest  to  test  whether  the  Rayleigh  parameters  (rms 

peak-to-trough  or  /~E)  of  the  20-rainute  stress  records — or  the  maximum  single 

stress  reversals  from  each  record  (X  ) — provide  a better  fit  to  a Weibull 

max 

distribution.  Nordenstrom  (16)  made  such  an  analysis  of  the  Wolverine  State 

data,  as  shown  in  Figure  10,  where  the  Weibull  distribution  of  both  / E and 

X is  seen  to  be  good.  Since  this  information  was  also  available  from 
max 

the  Esso  Malaysia  data,  plots  were  made  on  Weibull  paper  as  shown  in  Figure 
11.  As  in  the  case  of  the  stress  reversals,  these  data  plot  as  curves,  although 
the  Weibull  fit  is  better  for  /~E  than  for  the  stress-reversal  data  of  the 
Esso  Malaysia.  Similar  rms  peak-to-trough  stress  data  have  also  been  determined 
for  the  high-speed  SL-7  containership  Sea-Land  McLean,  and  the  distribution 
plotted  in  Figure  12,  indicating  a similar  nonlinear  trend.  (As  previously 
noted,  stress  reversal  data  were  not  available  for  this  ship). 

It  appears,  in  summary,  that  the  following  observations  can  be  made 
about  the  Weibull  distribution  of  ship  stress  data: 

1.  Weibull  fits  to  stress  reversal  data  are  poor . 

2.  Weibull  fits  to  / E values  are  fair  to  good . 

3.  Weibull  fits  to  X values  are  poor, 

max  


r 


Previous  conclusions  ( 17)  that  the  Weibull  distribution  of  / E values  gave 
excellent  agreement  with  statistical  data  were  based  largely  on  an  analysis 
of  Wolverine  State  data,  which,  as  we  have  shown,  seems  gratuitously  to  be 
very  well  fitted  by  Weibull  distributions.  When  applied  to  other  vessels, 
although  the  agreement  is  not  as  good,  it  appears  to  be  acceptable.  Previous 
work  (11)  shows  that  the  normal  distribution  of  /e”  values  is  also  acceptable. 

Grouped  Data,  Normal  JyT  Distribution 

The  method  of  extrapolating  ship  stress  data  to  predict  long-term  trends, 
which  has  been  the  standard  procedure  used  by  researchers  at  Webb  Institute, 
was  described  in  Part  1 of  this  report.  Briefly  repeated,  the  procedure  assumes: 

- Normal  distribution  of  / E values  within  each  Beaufort  number 

group,  with  observed  standard  deviation  for  each. 

- Rayleigh  distribution  of  peak-to-trough  stresses  within  each 

record. 

- Actual  distribution  of  Beaufort  Numbers  as  experienced  in  service. 

Then  three  integrations  are  performed  to  obtain  the  long-term  curve.  This 
procedure  does  not  require  any  specific  formulation  of  a distribution  function 
for  the  long-term  curve. 

For  the  four  ships  for  which  stress  histograms  were  available,  long-term 
distribution  curves  computed  in  this  way  have  been  published  previously  for 
the  Esso  Malaysia,  R.G.  Follis  and  Fotini  L in  (15 ) . The  results  are  shown 
in  Figure  13,  in  which  each  long-term  distribution  curve  is  plotted  along  with 
the  histogram  data  points  for  each  ship. 

A study  of  Figure  13  shows  that  the  fit  of  these  "grouped  Normal  / E " 
long-term  curves  to  the  actual  measured  data  is  good  for  the  R.G.  Follis 
and  excellent  for  the  other  two  vessels.  The  agreement  between  curves  and 
data  over  the  full  range  of  the  data  establishes  confidence  in  the  curve 
as  an  extrapolator  to  lower  probabilities.  That  is,  it  is  reasonable  to 
assume  that  the  probability  model  used  here  should  correctly  describe  the 
stresses  over  a long  period  (20  years  or  more)  since  it  correctly  describes 
the  cumulative  stress  distribution  observed  over  a period  of  two  or  three 
years,  provided  that  the  conditions  of  operation  (route,  weather  distribution, 
and  ship  characteristics)  remain  unchanged. 
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4. 
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The  Wolverine  State  analysis  was  published  in  references  (17),  (11),  (13 
and  (18) . A number  of  different  long-term  curves  were  calculated  from  these 
reports,  based  on  varying  amounts  of  data  available  over  the  years,  and  on 
modified  instrumentation  and  analysis  procedures  described  in  (18).  The  long- 
L term  curve  representing  the  most  comprehensive  data  population  is  from  (18), 

reproduced  here  as  Figure  14.  Two  long-term  curves  from  grouped  data,  assuming 
a normal  distribution  of  / E within  each  weather  group,  can  be  seen  to  fit  the 
histogram  data  well.  The  curves  are  based  on  20  and  44  voyages  respectively, 
while  the  histogram  represents  30  voyages.  The  Weibull  distribution  of  /""I 
(not  grouped)  shown  in  Figure  10  is  also  plotted  in  Figure  14,  and  the  fit  is 
also  acceptable.  Finally,  the  exponential  fit  to  stress  reversal  data  (see 
Part  1)  has  also  been  calculated  for  the  Wolverine  State.  This  distribution 
plots  as  a straight  line  on  Figure  14,  and  since  it  lies  everywhere  below  the 
measured  data,  it  is  considered  unacceptable  as  an  extrapolator . 

Of  the  two  methods  which  give  acceptable  fits  to  the  histogram  data, 
which  is  a better  or  more  reliable  extrapolator  to  much  longer  times?  The 
answer  to  this  question  is  somewhat  speculative,  because  verification  would 
require  a much  longer  data  collection  effort  to  extend  the  histogram,  and  the 
data  collected  on  the  Wolverine  State  represents  the  most  extensive  program 
of  its  kind.  In  principle,  grouping  the  stress  data  by  weather  severity  and 
analyzing  the  several  groups  as  independent  populations  should  be  more 
reliable,  because  the  scatter  of  /IT  values  within  any  one  weather  group 
would  be  less  than  that  of  the  entire  population  taken  as  a single  group. 

On  this  basis,  the  "grouped  normal  /e~"  technique  might  be  chosen.  For 
numerical  analysis,  the  normal  distribution  model  also  has  the  advantage  of 
having  parameters  (mean  and  standard  deviation)  which  are  easily  calculated 
from  the  data  without  necessitating  special  plotting,  as  is  necessary  in 
determining  the  Weibull  parameters.  Of  course,  the  weather  grouping  procedure 
can  be  combined  with  the  Weibull  distribution  assumption  by  describing 
separate  Weibull  fits  to  each  group  of  data,  as  is  done  by  Nordenstr^m  (19). 
This  procedure  would  have  the  advantage  mentioned  above  for  grouped  data, 
but  would  still  require  the  special  plotting  technique. 
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Conclusions 


In  the  first  part  of  this  report  it  has  been  shown  that  the  Weibull 
distribution  appears  to  be  a particularly  useful  tool  for  statistical  studies 
of  stochastic  processes.  Accordingly,  characteristics  of  the  distribution 
were  presented  for  convenient  reference,  its  estimation  discussed,  and  com- 
parisons made  with  a normal  distribution. 

It  was  shown  that  the  Rayleigh,  Voznesensky  and  exponential  distributions 


are  special  cases  of  a 

two-parameter 

Weibull 

distribution,  with  parameters: 

Rayleigh 

a = 1/R 

6 - 

2 

Voznesensky 

a = 1 

B = 

n 

Exponential 

a = 1/n 

B = 

1 

Finally,  the  generalized  Rayleigh  distributed  random  variable  is  in  reality 
the  sum  of  two  independent  random  variables,  one  of  which  is  normally  distributed 
while  the  other  is  Rayleigh.  As  the  band  width  varies,  the  generalized  Rayleigh 
distribution  tends  to  one  of  these  cases.  A 3-parameter  Weibull  distribution  may 
be  closely  fit  to  the  generalized  Rayleigh  distribution  with  the  Weibull  parameter 
dependent  on  the  band  width. 

In  Part  2 the  application  of  the  Weibull  distribution  to  the  problem  of 
ship  hull  stress  statistics  led  to  the  following  conclusions  regarding  the 
goodness  of  fit: 

Stress  reversal  data  (X)  - poor 

Rayleigh  parameters  (R  or  / E ) - acceptable 

Highest  stresses  per  record  (Xmax)  ~ poor 

Accordingly,  it  appears  that  the  numerical  solution  of  the  long-term 
distribution  is  most  appropriate,  using  either  a Weibull  or  normal  distribution 
of  Rayleigh  parameters.  It  is  believed  to  be  preferable  in  extrapolating  stress 
data  to  long  periods  of  time  to  classify  and  analyze  stress  data  by  weather 
groups. 
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Appendix  1 
The  Error  Function 

The  error  function  erf  (x)  which  appeared  throughout  the  main  body  of  this 
note  is  defined  as 


erf (x) 


= -2_  f o-'2 


dt 


The  full  list  of  properties  of  this  function  may  be  found  along  with  tables  in 
Abramowitz  (20).  The  most  commonly  used  properties  are  listed  below. 


lim  erf(x)  * 1 

X-*x> 


erf(-x)  - - erf (x) 


d erf (x)  _ _2  e *2 

dx  ft 

The  error  function  may  be  approximated  by  the  series  expansion 

erf  (x)  - ^ I 

/n  n-0  n!  <2n+1) 

The  error  function  is  also  associated  with  the  complementary  error  func- 
tion erfc(x). 

erfc(x)  s 1 - erf(x) 


2 

7* 


e~t2  dt 


we  also  have 


d erfc(x) 
dx 


- _2 

R 


_-x2 


erfc(-x)  « 1 - erf(-x)  - 1 + erf(x) 


Appendix  2 


The  Gamma  Function 


The  gamma  function  r(x)  is  defined  as 


/ tx-l  e-t 


This  integral  converges  for  x > 0,  and  via  integration  by  parts  satisfies  the 
recursive  relation 


p 00  - (X  - l)  r(x  - l) 

r(l)  ■ J e-t  dt  = 1, 


Since  r(l)  ■ J e-t  dt  = 1,  for  any  integer  m 

0 

F(m)  ■ (m  - 1)  ! 


T(l/2) 


1 t'1/2e-t 
0 . 

2 j e'u2  dl 


/n  erf  (®)  - /S' 


Additional  forms  of  the  gamma  function  are 


(t  - u2) 


r(x)  - lim  OJ — s!! 

n^ca  x(x+l)  . . . . (x+n) 


“ x -x/n 

x eYx  n [ (i  + f ) e ] 

n-1 


for  x ^ 0,  -1,  -2, . . . , 


where 


y m .5772156649. 


is  Euler's  constant. 


Although  the  form  of  the  gamma  function  is  rather  involved,  extensive  tables 
are  available,  (Abramowitz  (20)),  Quick  estimations  of  r (x)  may,  however,  be 
made  through  the  recursive  property  and  the  minimax  polynomial  approximation. 

T(1  + x)  - 1 - .577191652  x 

+ .988205891  x2 


C '' 


- .897056937  x3 
+ .918206857  xA 

- .756704078  x5 
+ .482199394  x6 

- .193527818  x7 
+ .035868343  x8 

which  represents  r (1  + x)  for  0 < x ^ 1 within  an  error  of  3 x 10  7. 

Extensive  additional  properties  of  the  gamma  function  may  be  found  in 
Abramowitz  (20) , Bateman  (21) . 
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Appendix  3 


Calculation  of  Moment  Generating  Function  for  the 
Generalized  Rayleigh  Distribution 


A random  variable  X is  said  to  have  the  generalized  Rayleigh  distribution 
if  its  density  function  is  x/l-t2 


f(x)  - — — [ e e~x2/2e2  + /l-  2 X /2 

/2n~ 


2 -x  f 
x e J 


e"t2/2  dt  j 


_o»  < x < 


The  change  of  variables  x » 2 u simplify  the  calculations,  and  if 

/1-e 

Mu(0)  is  the  moment  generating  function  of  u,  then  ^(0),  the  moment  generating 
function  of  x,  will  be 


Mx(0)  = M.eu  (0)  = m US*  ) 

vl-er 


The  probability  density  function  of  u,  g(u),is  f(x(u))  — *»£ 

V/2(l-£2)  . » 


or 


where 


and 


8<u)  ' e2  { 


ue 


-u2/2[e2/l-e2] 


Air  /l-e2 


e-c  /2  dt  } 


Mu(6) 


e®u  g(u)  du 


- e2[G(0)  + 


A^ 


H(e)  ] 


0(6)  - [ -"6  e-“2/2(l-t2) 

-w  Air  /1-e2 


du 


H(0) 


\J\ 


. u2  r _e2  , . t2 

u0  7 i-e2  ^ 


e ue 


dt  du 


The  integral  in  G(6)  is  quickly  recognized  as  the  integral  defining  the 
moment  generating  function  of  a normally  distributed  random  variable  with 


Thus 


p = 0 and  o2  - 1 - e^ 
G(o)  - e 2 
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The  integral  for  H(6)  is  more  involved.  If  the  order  of  integration  is  reversed 


H(6) 


-1 


0u 

ue  e 


2^  1 T-e2  J 9 

e du  dt 


which  is  of  the  form 


w U 

JI 


-l/2l-^  + -j|] 

oiz  022 


ye®^  e 


dy  dx 


V o * w 

where  y = u,  x = t,  B 1,  t»2  = ”^2 

By  completing  the  square  in  the  integrand  we  obtain 


1-e2 


I - e 


52o22 


v - 0o22  -2 
3 „ 2 ] 


dy  dx 


- 00  X 


or 


■ ■ J 1 


[ y + 0OO  ] e 


-1/2  [ ^ ^ ] 

o2  o22 


dy  dx 


x-002^ 
2 

where  y has  replaced  y - 002 
Thus  I « Ij  + I2  where 

-y  00  00 

11 


-1/2  [ + 2^  ] 

0 022  e °1  °2  dy  dx 


x-0q 


'2 


and 


I2  - e 


e2oz2 


/oo  00  9 

’•  1 I . .-“"i? 


+ ^ 1 


9 

-00  X-0O2 


°22  dy  dx 


We  now  examine  1^  and  I2  and  evaluate  them  through  several  changes  of 
variables.  Dnder  a change  of  variables  x = f(t,s),  y = g(t,s)  the  integral 


J | F (x,y)  dx  dy  = JJ|  F(f(t,s),  g(t,s))  | 


dt  ds 


whenever 


3(t,s) 


* 0 


where  R*  is  R expressed  in  terms  of  the  new  variables  and 


3x 

d (x.y)  _ 

d(t,s) 

3y 

ds 

ix 

3t 

ds 

First  let  t « ~r  , s * 2—  • 

01  o2  3(s,t)  °1°2 


and  becomes 


e2°272 


h = 6 


e a,o. 


-t2/2  -s2/2 


ds  dt 


“ a, t 

x 02 


We  now  relate  the  t,s  coordinate  system  so  as  to  make  the  line 
o 1 1 

s = - 002  parallel  to  the  horizontal  axis.  Thus 

0o 


t • t ' cos  0 - s ' sin  0 
s = t1  sin  0 - s'  cos  0 


where  tan  0 = — , sin  0 = 

o2 

Hence 


/o^2+022 


, cos  0 


'a  2+o  2 
1 2 


[o2t  ’ - OjS  »],  s 


[Ojt  ’ + o2s  ']  , 


Jl  °2 


3 (s 1 , t ' ) 


1 and  e 


t2+s2 

2 


:')2+  (s')2 


If  we  substitute  the  primed  variables  into  and  the  rewrite  the  integral  with 
the  primed  notation  eliminated  we  have 


02O22/2 


00  _«-2 


3 r -t2/2  r -SV2 

!02  t J 6 dt]  [ J e ds  ] 


-s2/2 


7 -t2/2 

I 6 dt  ■ 

-00 

“ -s2/2 

e ds 

0 o22 

/a^  + 022 


• 


f -u 

e du 

' e a22 

/l  i/a^2  +o~^ 

f° 

1/2  J e U du  + *^2 

-00  0 


6 o2* 


/2  /a^2  + o2 2 


— u2  , 
e du 


Hence 


. ^2  j£_  erf  ( 

2 2 


9 o92 

2 ) 


Jl  /a22  + o22 


/f  j^+  erf 


/2  ^ 


2 

==r) 
+ C22 


°WI 2 


6 o2~ 


•E 


+ erf  (- 


9 02* 


/?  /aj2  + 02Z 


Next  we  consider  I2  under  the  same  set  of  variable  changes. 


I2  = e 


i '■ 


92o22/'2 


-1/2  [ + -4  1 

ye  °1  CT2  dy  dx 


—00  9 

x-0<72z 


= e 


t 


9 2o22/  2 


a a 


s 0 ds  dt 


fil  - 9o2 

a2 


°1  a2 


^7  + o2 


2 ee2°22/2  r -] 

■ ■ --  U2I  + 122J 


where 


[21  " °1 


-1/2  [t2+s2] 
t e ds  dt 


” 9 a22 


°1 


*^ai2  + 022 


f -t2/2 


t e dt 


S’  \ 

OO 

-1/2  s2 


u 


ds 


9 02"“ 


/al2  + a22 
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- 


0 


since 


f -t2/2 

t e dt 


while 


But 


22  = °2 


IJ 


-1/2  [t2  + s2] 
s e ds  dt 


—co  0 

0 02^ 


JZ 22  + o22 


r -t2/2 

f f -s2/2 

!22  - cj2 

11  * 

J J s e ds 

1 

o22 

/o^2  + o22 


2 4 

*•  rt  r\  ~ 


0Z  02 


/27 


°2 


2(o12  + o22) 


Thus 


02o „2/ 2 


2/2 


01  o2' 


/o^2  + a22 


e2  o2" 


0 + o2  e 2(o2"  + ^ J 


Finally 


I = e 


02a22 


/2 


6 ox  o23  Tf  J 1 + erf  ( 


0 o2* 


°1  q2 
/o^2+  o22 


02a22 


2/oi2+o22 

2 / q22 

j 2 o 


1 


If 


°!2  " 1 » a22* 


then  a j2  + o22  = 1/e2 


and 


e2/2  ( icp  J 

H(0)  - e e 


0Or^-3/2  _ , [ ! + erf  ( e_  Jbd  ) ] 


+ /27  e 


3/2 


02 [1-e2]  1 


h 


/2  *= 
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FIGURE  1 R(B)  vs. 


Fig.  3— b : Normal  and  Weibull  Distributions  of  Ship  Bending  Stress  Data,  Weather  Group  III 


Fig*  3-c:  Normal  and  Weibull  Distributions  of  Ship  Bending  Stress  Data,  Weather  Group  IV. 


Fig.  3-d:  Normal  and  Weibull  Distributions  of  Ship  Bending  Stress  Data,  Weather  Group  V. 


o.qssq'? 


Jo.'ns’* 
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Fig.  7 : 


i.o.  2.0  3.0  H-0  So  to  T.o  S.o  TJ)  i&o 

5Tf?eSS  (HPSl) 

Esso  Malaysia  Stress  Reversal  Histogram  on  Weibull  Probability  Paper, 
with  Possible  Weibull  Fits.  » 


• • I I 


1.0.  2.0  3.0  H.o  So  to  T.O  f.O  IjO  |flO 

STRESS  f KFS\ 

Fig.  9 ; Stress  Reversal  Histograms  on  Weibull  Probability  Paper 
from  Measured  Stresses  on  Several  Ships 


Fig.  10:  Weibull  Predictions  of  J E and  Xmax  Distributions, 
Wolverine  State  (From  it) 


Fig.  11:  Distributions  of  /e"  and  Xmax,  Esso  Malaysia 
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Fig-  12'  Distribution  of  /lT,  Soa-Land  McLean 
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In  Part  2,  the  most  promising  long-term  distributions  derived  from 
the  above  survey  are  applied  to  ship  bending  stress  data  from  four  ships. 
It  is  concluded  that  a numerical  solution  of  long-term  distributions, 
using  either  Weibull  or  Normal  distributions  of  the  short-term  Rayleigh 
parameters  (classified  by  weather  groups)  is  better  than  any  explicit 
function.  ^ 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfUTian  Data  Enl.rad; 


